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1. Introduction 

A major focus of interest in hadronic physics is the quest to understand from first principles 
the structure of mesons and baryons. In particular form factors yield information about the size and 
shape of the hadrons. 

Much theoretical and experimental work has gone into understanding the structure of nucleons 
and the N — > A transition. Lattice QCD calculations of the nucleon and N — > A form factors (FFs) 
have been carried out within the same lattice setup as the one used in this work ^, |p. Exper- 
imental information on the FFs of A(1232) is scarce due to its short lifetime (~ 10~ 23 s) ^J. 
However, in a finite-volume simulation with heavy pions, the A is stable and accessible to lattice 
techniques. A pioneering lattice study [Q] investigated the electromagnetic form-factors of the A in 
the quenched approximation. Recently, a state-of-the-art lattice calculation of the electromagnetic 
FFs of the A and the associated transverse charge densities in the infinite momentum frame has 
been carried out rf7h. In this report we extend the effort to the axial and pseudoscalar form factors 
of the A(1232) and present preliminary results. To our knowledge this is the first time that these 
FFs have been computed. 

Despite the difficulty of experimental confirmation of these results, they can yield an evalua- 
tion of the axial charge and the effective 7iAA couplings, parameters which can be fed into chiral 
expansions to aid the chiral extrapolations of, for example, the nucleon axial charge. The axial 
Ward-Takahashi identity relates the axial FFs to the pseudoscalar FFs. As in the nucleon case, one 
can derive the generalized Goldberger-Treiman relations. In this work we derive these and check 
their validity. 

2. Lattice calculation 

This project closely follows the methods used for extracting A + electromagnetic form factors 
as described comprehensively in Ref. [R]. We begin with the expression of the isovector axial 
vertex: 

(A(p f ,s f )\A^\A(p h s i )) = u a ( Pf ,s f ) [^ A ] aP u p ( Pi , Si ), (2.1) 

with 

A^(x)=W(^Y5yV( x )- (2-2) 
The right-hand side is an expression containing the most general decomposition of the axial 
vertex in terms of four form-factors, which we label g\, g%, h\ and hy. 

^ A = - 8 ap (^V^+rf^y 5 ) +^ (M^yV+^lj-T 5 ) , (2-3) 

and u a is the Rarita-Schwinger spinor and q = pf — pt. 

Similarly we can write the pseudoscalar vertex in terms of two form-factors, g and h: 

(A( Pf ,s f )\P\A( Pi ,s i )) = u a (p f ,s f ) [0 ps ] ap u p ( Pi , Si ) (2.4) 

with 

P(x) = \if(x)r 5 — Y(x) (2.5) 
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and 



4Ma 



(2.6) 



We isolate the form-factors by constructing ratios of lattice two- and three-point functions. The 
standard lattice interpolating field with the A + quantum numbers is given by 

Xt + a (x) = -Le^ [2 (u aT (x)C 7(y d b (x)) u c a (x) + (u aT (x)C 7a u h (x)) d c a (x)] . (2.7) 

The two-point and three-point functions of interest are: 

G X ^(T v ,q,t) = Z e ^^aiXaa(tf,Xf)X^t,m,A0M 

X,Xf 

G a ,(r v ,p,t) = Ze-^Tl, a (Xaa(t,x f )Xra^M 



(2.8) 
(2.9) 



where X stands for the axial- vector and pseudo-scalar currents defined in Eqs. (2.3) and (2.5) 
respectively and 



1 



r 4 = -(i + /), r* = -(i + /))$%, * = 1,2,3. (2.10) 

We examine ratios of these to eliminate unknown Z-factors and leading time-dependence.: 



G x a ^(T,q,t) / G kk {T\p h t f -t)G kk {Y\0,t)G kk {T\0,t f ) 



G kk (r 4 ,0,t f ) V G kk {Y\0,t f -t)G kk {T\p h t)G kk (T\p h t f ) 
These ratios tend to a constant at large Euclidean time separations tf — ti and t: 



(2.11) 



-ftff(/i)T(r,^,0 z f — 

with the kinematical constant given by 



cn 



Ctr 



C 



M A Ml 



Ml 



2E A( Pi ) 2E A( Pl ) E A( Pl ) , E A(p.) 



Ml 



(2.12) 



(2.13) 



3. Lattice simulation 

We utilize 200 quenched Wilson configurations on a lattice of size 32 3 at j3 = 6.0, which cor- 
responds to inverse lattice spacing of a -1 = 2.14(6) GeV. We perform the analysis at three hopping 
parameter values K = 0.1554, 0.1558 and 0.1562, corresponding to pion masses m n = 563,490, 
and 411 MeV respectively. Additionally, we use a mixed action of domain wall valence fermions 
on a staggered sea simulated by the MILC collaboration with an Asqtad improved action []|], with 
a pion mass of 353 MeV. A total of 200 configurations are analyzed at one value of the pion mass. 
The details of the simulations are summarized in Table [j] 

We use the sequential source method |Q] to calculate three-point functions. We use the same 
fixed source-sink separation as was used in Ref. ^ of ~ 1 fm, or 11 time-slices on the three 
quenched ensembles and eight time-slices on the mixed-action ensemble. 
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Wilson fermions 


V 


# confs 


K m n m n /m p mpi 
(GeV) (GeV) 


(GeV) 




SIM-I 


: Quenched, J3 = 6.0, a~ l = 2.14(6) GeV 




32 3 x 64 


200 


0.1554 0.563(4) 0.645(9) 1.267(11) 


1.470(15) 


32 3 x 64 


200 


0.1558 0.490(4) 0.587(12) 1.190(13) 


1.425(16) 


32 3 x 64 


200 


0.1562 0.411(4) 0.503(23) 1.109(13) 


1.382(19) 






SIM-II: Mixed action 




Asqtad ( 


am u,dls 


0.01/0.05), DWF (am u4 = 0.0138), a' 1 = 1.58(3) GeV 


28 3 x 64 


200 


0.353(2) 0.368(8) 1.191(19) 


1.533(27) 



Table 1: Ensembles and parameters used in the calculation of form factors. 

4. Extracting form-factors 



For each ratio 



IT 



o(H)t 



tr 



rA 



GO" 



we work out the trace algebraically for specific combinations of a, T, and T ; or r 4 : 

3 3 

V' '■'<" Z- U i o^[T J A (ja i{pf)ff a/ ^A T ^{pi)) 



or 



with 



3 

I to 

C\T=1 



•tr [^haa'iPfWa'^'^ipi)} 



1 
1 
1 



and 





-1 
1 



1 


-1 



(4.1) 



(4.2) 



(4.3) 



(4.4) 



We refer to these traces as Type I and Type II respectively. For the right-hand sides we now have 
linear combinations of the form-factors where the coefficients are functions of the initial energy, 
Ma and the spatial initial momentum (We Wick rotate and work in the rest-frame of the sink). 
In general, the form of the expression is different for jU = 4 and /I = 1,2,3. The left hand side is 
calculated on the lattice, so we may now solve a system of linear equations to isolate the form- 
factors. The FFs are extracted by the simultaneous over-constrained analysis of all the relevant 
ratios that contribute to the transition per given Q 2 . The renormalization constant Za is required for 
the axial FFs. They are given in Table IV in [Bp. 

We summarize the results obtained for the axial form factors calculated on all four ensembles 
in Figure [l[ By extrapolating the g\ curve to Q 2 = we get an estimate of the axial charge of the 



A + . In order to connect gi (0) to the axial charge #aa> we use the relations given in [10]. We find 
g A = 1.2694(28) [flip we expect that 



that g\ (0) = — |gAA- Using the tree-level SU (4) relation gAA = — fg/t and the experimental value 



ft (0)«~(lJ7) 



0.76, 



(4.5) 
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Figure 1: Results for the four axial form-factors, g\, g3, hi and /13. The dashed curve shows consistency of 
gi for the K = 0.1562 ensemble with a dipole fit giving a pole mass of 1.67(5)GeV. 



which is consistent with our results. 



5. Effective couplings 



Referring to Eq. (|2.4|), we decompose the A matrix elements of the pseudoscalar current into 
two axial 7iAA couplings, G x aa and H k aa, with the relation: 



2m q (A Pf \P\A Pi 



and we identify: 



\E A (p f )E A (pi) J (ml-q 2 ) 



g afS G nAA (q 2 



> + 



H 



7TAAI 



t "fua (5.1) 



m q g 



fnmlGxAAiq 2 



and 



(mt 



■r 



m q h : 



f n m 2 H nAA (q 2 



(mt 



■r 



(5.2) 



Note that because the two different possible contractions of the Dirac indices of the 3/2-spinors 
give us two pseudoscalar form-factors we get two effective axial couplings, unlike the case of the 
nucleon and the N — A transition. The quark mass m q is computed from the axial Ward-Takahashi 
identity and f n from the pion-to-vacuum amplitude. Both are computed from appropriate combi- 
nations of two-point functions as shown in the reference P]. The renormalization factor Zp is not 
required as its occurrences in m q and g or h cancel. 

The results for these form-factors are shown in Fig. |2[ As can be seen, despite the large 
statistical errors, G k aa increases with decreasing Q 2 for the unquenched lattices. 
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Figure 2: The effective couplings G^aa an d H kAA for the quenched ensemble with k = 0. 1554. 
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Figure 3: Ratios of axial form-factors, g3/gi and h^/h\. Dashed curves show consistency with a pion pole 
fit for the K = 0.1562 ensemble. We omit the mixed ensemble from the hj,/h\ plot as the signal is washed 
out by the larger error bars. 



6. Goldberger-Treiman relations 

From the axial Ward-Takahashi identity, we have the relationship 

(A|^|A) = 2m,(A|P|A). 



Applying the momentum operator on (gjj) and ( |2.3| ), the left-hand side gives 

i a n?> 



(A|d M A"|A)=2m A 



l -Tg3)g aP + (h-Th 3 ) 



Ami 



UaTU 



with T = (27^- Usin § Ec l s - <0 we § et 



2m A ( 



2f Jt mlG 7tAA (q 2 



and 2mA (^1 — T/13) 



2/ 3r «4// 3rAA ( (? 2 



(6.1) 



(6.2) 



(6.3) 



If we demand that the #3 and terms cancel the pole at q 2 = m 2 , we get the Goldberger- 
Treiman relations. In Figure [3] we show that the ratios g$/g\ and h^/hx are consistent with pion- 
pole behavior. We note that, as with the effective axial couplings, there are two Goldberger-Treiman 
relations, namely: 



fnG n ^(q 2 ) = m A gi( 



and 



7tAA\ 



m/Ji\ 



(6.4) 
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Figure 4: The two Goldberger-Treiman relations for the A baryon. 



In Fig. ^| we plot the ratio of the left-hand to right-hand sides of the expressions given in 
For low Q 2 , these quenched ratios deviate from unity but are in agreement with unity for 
Q 2 ~ 0.8GeV 2 . Similar behavior was observed for G^nn and G^na in 01- In the unquenched 
ensemble the errors are too large to enable any conclusions. 

7. Conclusions 

In this work we have evaluated, for the first time, the A + axial form factors, g\, g%, h\, hi, 
as well as the pseudoscalar form factors g, h. We have shown that these axial and pseudoscalar 
vertex compositions yield two effective 7iAA couplings, G^aa, and H^aa, which in turn satisfy dual 
Goldberger-Treiman relations. Results obtained in the quenched theory are accurate enough to 
enable a check of these relations and show that there are deviations for small Q 2 values where chiral 
effects are expected to be large. Unquenched results using a mixed action have large statistical 
errors and require further analysis for allowing a definite conclusion to be reached. 
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